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Abstract. The observed variability of BL Lac objects 
and Quasars on timescales 1 day (intraday variability, 
IDV) have revealed radio brightness temperatures up to 
Tb ~ 10 16 — 10 20 K. These values challenge the beaming 
model with isotropic comoving radio emission beyond its 
limits, requiring bulk relativistic motion with Lorentz fac- 
tors r <; 100. We argue in favor of a model where an 
anisotropic distribution of relativistic electrons streams 
out along the field lines. When this relativistic beam is 
scattered in pitch angle and/or hits a magnetic field with 
components perpendicular to the beam velocity it starts 
to emit synchrotron radiation and redistribute in momen- 
tum space. The propagation of relativistic electrons with 
Lorentz factor 70 ~ 10 2 — 10 4 reduces the intrinsic vari- 
ability timescale At' to the observed value At ~ At'/ 70 
so that the intrinsic brightness temperature is reduced 
by a factor of order ~ l/7o, easily below the Inverse 
Compton limit of T h <, 10 12 K. When looking at a single 
event we expect the variability time scales At to be inde- 
pendent of frequency for a monoenergetic electron beam, 
whereas for a beam with a spread out distribution of en- 
ergies (e.g. power-law) parallel to the magnetic field the 
timescales are shortened towards higher frequencies ac- 
cording to At oc i/~ . The observations seem to favor 
monoenergetic relativistic electrons which explain several 
properties of variable blazar spectra. The production of 
variable X- and gamma-ray flux is briefly discussed. 

Key words: Acceleration of particles - Plasmas - Ra- 
diation mechanisms: non-thermal - BL Lacertae objects: 
general - Galaxies: jets - Quasars: general 



1. Introduction 

It has become a reliable observational fact that active 
galactic nuclei are variable on all timescales throughout 
the whole electromagnetic spectrum from the radio range 
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up to the gamma-ray-range (see Wagner and Witzel 1995 
for a comprehensive review). 

An especially interesting type of variability arises in 
the radio range on the timescales of less than one day, 
called intraday variability (IDV) . IDV at radio frequencies 
has brought up significant difficulties with the standard 
picture of beamed intrinsically isotropic incoherent syn- 
chrotron emission. If the variability timescale At is used 
to calculate the source size by I ~ cAt, where c is the ve- 
locity of light the resulting brightness temperatures reach 
up to the range Tb = 10 16 — 10 20 K, far in excess of the in- 
verse Compton limit of Tb = 10 12 K. The question arises: 
What is the physical nature of this very intense variable 
radio emission? 

We first summarize the properties of the varying 
sources. They usually vary on the 5-10% level at cm wave- 
lengths (Witzel 1992) but it may also be as high as 35% 
(Quirrenbach et al. 1992). The sources that show IDV of 
the total intensity generally also show IDV of the polarised 
flux and the polarisation angle, usually (anti-) correlated. 
The typical Blazar that is variable on such short timescales 
seems to have the following properties (Wegner 1994): 1.) 
it reveals a compact structure on VLBI scales (VLA com- 
pactness is not sufficient), 2.) associated variability of the 
polarisation, 3.) if the object is variable at radio frequen- 
cies, it is also variable in the optical, 4.) all BL Lac ob- 
jects are variable, 5.) sources with a positive spectral index 
a® ~ 0.2 between 2 and 6 cm (/„ oc v a ) show the largest 
variability amplitudes. 

In the standard picture of beamed intrinsically 
isotropic emission it is found that a brightness temper- 
ature T' in the comoving frame will be seen as a van- 
ability brightness temperature of Tb ~ D 3 T' , where 
D = [T(l — /3cos6>)] -1 is the Doppler factor of the jet 
(Blandford et al. 1990). Here T = (1 - /3 2 )" 1 / 2 is the bulk 
Lorentz factor, (3c the jet velocity in the observers frame 
and 9 the angle between the jet and the line of sight. The 
Doppler factors calculated from the observed superlumi- 
nal expansion arc of the order of D 10 (Cawthorne 
1991). This is much too low to explain the extreme val- 
ues of Tb found in sources with IDV which would require 
D <; 100. Begelmann et al. (1994) discuss the energetic 
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constraints on highly relativistic jets (T ~ 30 — 100 to ex- 
plain 7b IS 10 16 K) and find that the jet must carry kinetic 
energy fluxes comparable to the total power of a luminous 
quasar (L ^ 10 47 ergs -1 ), whereas many objects with IDV 
are low luminosity AGN of the BL Lac type. 

A completely different point of view involves coher- 
ent emission (Baker et al. 1988; Krishan and Wiita 1990; 
Benford 1984, 1992; Lesch and Pohl 1992). Coherent ra- 
diation easily accounts for huge Xb, since the radiation in- 
tensity is due to collective emission of all particles within 
a coherent volume. The rather involved physics is well 
known from plasma laboratory experiments. The central 
argument against coherent radiation processes in AGN is 
the limitation of the brightness temperature by induced 
Compton scattering and/or Raman scattering (Coppi et 
al. 1993; Levinson and Blandford 1995). However, it was 
shown by Benford and Lesch (1998) that the saturation 
arguments are not valid for strong Langmuir turbulence 
excited by relativistic electron beams and radiated coher- 
ently via inverse Compton scattering. 

In this contribution we would like to present an al- 
ternative scenario, which considers the possibility that 
the energy distribution function of the synchrotron emit- 
ting relativistic electrons is not isotropic in pitch angle 
throughout its propagation in the very center of an AGN. 
In other words we consider an electron beam, whose mo- 
mentum parallel to the magnetic field is considerably 
higher than its momentum perpendicular to the field 
lines. The advantage of this model is that the synchrotron 
brightness temperature of a beam distribution has to be 
calculated with the square Lorentz factor of the particles 
(instead of the bulk Lorentz factor for an isotropic distri- 
bution), which can be of order 10 2 — 10 4 . Then the dis- 
crepancy of the inverse Compton limit of 10 12 K with the 
observed values up to 10 20 K is resolved. 

In the next section we consider the formation and 
propagation of relativistic electron beams within the en- 
vironment of active galactic nuclei. Here we especially 
investigate the conditions for isotropization by magneto- 
hydrodynamic turbulence and how this redistribution to- 
wards larger pitch angles can be avoided in regions with 
high magnetic field strength. What follows is a calcula- 
tion of the resulting synchrotron radiation of beams for 
two cases: for a beam with a narrow energy distribution 
(quasi-monoenergetic with a clearly defined low and high 
energy cutoff) and a beam with a broad energy distri- 
bution (power law). Finally we discuss our results and 
its implications for different parts of the electromagnetic 
spectrum from the radio up to the gamma-ray range. 

2. Formation and propagation of a relativistic elec- 
tron beam 

In this section we consider the formation of a beam of 
relativistic particles near the central supermassive black 



hole in an active galactic nucleus and its propagation in 
the jet. 

2.1. Acceleration of beam particles 

We consider a spinning back hole with mass Mbh = 
10 8 MqMs and angular velocity f^BH with an accretion 
disk that brings in ionized gas and magnetic fields with 
strength B ~ lO 4 ^ G near the hole. The rotating mag- 
netic field lines excert a torque and accelerate particles to 
high energies (Gangadhara and Lesch, 1997). The Lorentz 
factors obtained in this way arc limited by the inverse 
Compton interaction of the relativistic particles with the 
ambient radiation field, basically UV photons from the 
disk. This leads to the production of X- and 7-rays. The 
electrons will lose their energy perpendicular to the mag- 
netic field very rapidly due to synchrotron radiation in a 
time 

is = 5 10" 3 7 a3 1 B 4 r 2 sin" 2 a sec (1) 

where a is the angle between the particle velocity and 
the magnetic field and 70 = 10 3 7o,3 is the Lorentz fac- 
tor of the relativistic electrons. The combined action of 
acceleration and inverse Compton losses results in a mo- 
noenergetic distribution of relativistic electrons propagat- 
ing along the magnetic field. The particles accumulate at 
the energy where the timescales for gains and losses are 
equal and balance each other. 

Another possibility for the formation of relativistic 
particle beams along the magnetic field is by means of 
reconnecting magnetic fields in the inner part of the accre- 
tion disks and their coronal environments. Schoppcr et al. 
(1998) have clearly shown that acceleration by magnetic 
reconnection and energy losses via synchrotron radiation 
and/or inverse Compton scattering leads to beam- like en- 
ergy distributions. Their test-particle simulations use re- 
alistic magnctohydrodynamical magnetic field structures 
originating in sheared, magnetized gas flows as the back 
stage for the acceleration process. Keplerian shear of the 
field lines leads to the build up of strong magnetic gradi- 
ents, corresponding to high current densities, which finally 
are dissipated by localized three-dimensional magnetic re- 
connection structures with a length d. Three-dimensional 
magnetic reconnection is equivalent to the presence of 
magnetic field aligned electric fields E\\ along which par- 
ticles are efficiently accelerated on very small time scales 
(Schindlcr et al. 1991). Thus, the maximum energy eE»d 
is reached by almost all particles, only a few particles are 
influenced by energy losses. Finally a quasi-monoenergetic 
electron distribution function results. 

We note that in the central 50 pc and the inner 10 13 cm 
of the Milky Way several locations (Arc and Sgr A* ) re- 
veal inverted radio spectra, in which the flux increases 
with frequency with an exponent 1/3, which is easily in- 
terpreted as optically thin synchrotron radiation of quasi- 
monoenergetic electron distributions (Lesch and Reich 
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1992; Duschl and Lesch 1994). Such a spectrum has also 
been observed in M81, which is classified as a "dead" 
Scyfert nucleus (Reuter and Lesch 1996; Bottcher et al. 
1997). 

Obviously central regions of galactic nuclei are pre- 
ferred physical environments for the formation of beams. 
The main reason is that both mechanisms, centrifugal 
acceleration and magnetic reconncction are natural con- 
stituents of a rotating magnetized system, including tur- 
bulent gas motions and directed flows (jets), which are all 
supposed to be necessary ingredients of models for active 
galactic nuclei (e.g. Blandford ct al. 1990 for a general 
review) 

2. 2. Stability of the beam 

The question arises if these streaming particles can ex- 
cite Alfven waves and if they are scattered in pitch angle, 
which would lead to isotropisation. It is generally believed 
that this occurs throughout the whole jet. Instead we will 
argue now that this can not be the case for the inner- 
most part considered here. The phase velocity of Alfven 
waves with frequency w and wavenumber k propagating at 
an angle 9 with respect to the magnetic field is given by 
(Krall and Trivelpiece 1973) 



k 2 



Vl cos 2 6 



Here the parameter called the Alfven velocity is 
V A = B/ ^/AirmpUp 



(2) 



(3) 



where m p is the mass and n p is the number density of 
the protons constituting the background plasma that car- 
ries the waves. The numerical value is given by Va = 
2.2 10 11 n p X ^ 2 B cm/s. Near the black hole this is larger 
than the speed of light for n p < 5.3 10 9 £>f cm~ 3 . If one 
asssumes that the total kinetic power of the jet is in ther- 
mal material this is just about the maximum average gas 
density in the jet (Celotti et el., 1998). Since it is much 
more likely that the gas is confined in filamentary struc- 
tures (probably on the rim of the jet), the density of gas 
that can carry Alfven waves in the central regions of the 
jet is much lower. In fact from the lack of significant Fara- 
day rotation in BL Lac objects Celotti et al. (1998) de- 
rived rip/v 10 6 cm -3 at the base of the jet (with radius 
i?o = 10 14 cm) where fy is the volume filling factor of the 
thermal material. So either the density or the the volume 
filling factor is low. They concluded that most of the en- 
ergy in the jet is carried by the relativistic particles and 
the magnetic field in BL Lac objects. 

The relativistic electrons can generate and interact 
with the waves when the resonance condition 



uj — sCt — k\\v\\ 







(4) 



is fulfilled. Here k\\ = kcosO and v\\ — vcosa are the com- 
ponents of the wave vector and the particle velocity along 
the magnetic field, respectively, with a being the pitch 
angle. Further, s is the harmonic number of the gyromag- 
nctic interaction and Q = Q e /70i with Q c = eB/mc, is 
the relativistic gyrofrcquency of the electrons. The case of 
resonant scattering of Alfven waves occurs when (Melrose, 
1986) 



|fc||i>l|| « \sQ\ > u> . 



(5) 



From equations (4) and (5) it follows that resonance is 
only possible when 



1/2 



+ /?cosa>l 



(6) 



where (3 is the velocity of the particles normalized by c. In 
the case Va ^> c this condition reads (3 cos a 3> 1, which 
cannot be fulfilled. When the phase velocity of the Alfven 
waves along the field is close to c the waves that prop- 
agate in the medium are almost purely electromagnetic; 
these are ineffective in pitch angle scattering. Achatz and 
Schlickciser (1993) find that an electron-positron beam is 
completely stable for 



< 



To 



1/2 



(7) 



where lu p = (47rn c e 2 /m c ) 1 / 2 is the plasma frequency. 
This can be compared with the condition as discussed 
above, since we can write Va/c = (m c /m p ) 1 / 2 ri c /cijp <; 
(7om c /m p ) 1 / 2 ~ 1. So for the Lorentz factors considered 
here both are essentially the same. 

It seems that in the acceleration region near the black 
hole anisotropy is favored and a quasi monoenergetic and 
one-dimensional beam of relativistic particles can propa- 
gate along the magnetic field. 

2.3. Disruption of the beam 

The relativistic beam distribution of electrons with num- 
ber density N and 



N(j) = N (5( 7 - 70 ) 



(8) 



and the background plasma carry a current that produces 
a toroidal component B± of the magnetic field near the 
rim of the jet. The strength scales with the radius R of 
the jet as 



Bj_ = Bj_,o ( -gr 



(9) 



and confines it magnetically. The poloidal magnetic field 
in the central part of the jet falls off as 



B» =B 



||,o 



R_ 



(10) 
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As the jet widens in radius the parallel field ceases to be 
important. 

It is now interesting to ask when the Alfven velocity 
Va becomes less than c. The scaling (9) for the toroidal 
magnetic field gives a constant value of Va when the par- 
ticle density scales according to 



R 

Rq 



(11) 



But since the internally dominating poloidal field falls off 
more rapidly we find that the Alfven velocity is below the 

1 /2 

velocity of light for B/n p & 0.1 or 



R_ 

Ro 



Z 10 B 0A n 06 



1/2 



(12) 



with no = 10 6 no,6 cm~ 3 . This means that the beam of rel- 
ativistic electrons will not be scattered in pitch angle as 
long as the jet has not expanded enough. When condition 
(12) is fulfilled, i.e. when R <; 10 16 cm and consequently 
B± ~ 10 2 G the beam will rapidly start to excite Alfven 
waves which scatter the particles in pitch angle leading to 
isotropisation. The flow of the particles becomes unstable 
and turbulent as the relativistic beam particles couple to 
the background medium. A standing shock will develop 
some distance L away from the black hole, which can be 
identified with the VLBI core. The relativistic electrons 
then run into a shock compressed magnetic field with large 
perpendicular components and impulsively start to radi- 
ate. 

3. Emission from the relativistic electron beam 

The radiation will be dominated by the perpendicular 
magnetic field component. The relativistic particle beam 
hits the magnetic field with a Lorentz factor 70 and starts 
to emit synchrotron radiation. Although the individual 
electron will have its emission confined to an angle of the 
order of I/70, the overall emisson of the jet will be spread 
over an angle tp ~ 0.1, the opening angle of the jet. The 
spectrum reaches up to the frequency (e.g., Longair 1981) 



where At' is the corresponding time interval in the beam 
frame and V is the intrinsic length scale of the beam. In 
the case of ~10 hrs variability from a 70 = 10 3 beam 
the length scale is given by V ~ 10 18 cm. The brightness 
temperature that is calculated from light travel time ar- 
guments is then overestimated by a factor 7 2 since 



cx (cAt)- 2 = (cAt'/ lQ )- 2 = r 2 7o 2 



(15) 



and I' is the comoving size of the emitting region. The 
variability timescale is the same at every frequency, from 
the optical down to the radio. The observed variability 
timescale will be given by (14) as long as it is larger than 
the radiative loss time, see eq. (1). If the latter becomes 
larger it smears out the first, because then the whole beam 
seems to flare up with the synchrotron life time of the 
particles. From eqs. (1) and (9) we find the condition 



A^3io3 7 V 2 Bo4 



At 
10 hrs 



1/2 



(16) 



assuming B± to be dominant during the emission, accord- 
ing to what is seen in BL Lac objects. As long as the ra- 
dius is limited in such a way, variability is dominated by 
Dopplcr boosting. Indeed jet formation models yield fairly 
constant radii over large distances (Camenzind 1996). 

If we assume the energy distribution of the beam to be 
a power-law (still with zero perpendicular momentum) 



N(j) cx 7 



(17) 



the emission of synchrotron radiation at a certain fre- 
quency will be dominated by the particles with a Lorentz 
factor given by relation (13). The variability timescale 
then depends on the frequency. Using (13) and (14) the 
dependence is found to be 



At cx 7 cx v 



-1/2 



(18) 



This means that for a spread in the energy distribution of 
the radiating particles the spectral emission will be vari- 
able on a shorter timescale at higher frequencies. 



v c = ^To^e sin a = 4.2 10 7 3 B 2 sin a Hz 



(13) 



where v e = £! c /27r is the nonrelativistic gyrofrequency. 
Thus we have emission from radio frequencies up to the 
optical. The intensity I v is proportional to z/ 1 / 3 for fre- 
quencies v <C v c and falls off exponentially for v v c . 
This naturally gives the observed inverted spectra of com- 
pact radio sources in Blazars. Self- absorption effects may 
become important at lower radio-frequencies. 

The duration of a radiation event at a fixed location 
in the observers frame is given by 



At 



7o 



4. Further evolution of the spectrum 

When the beam starts to dissipate its energy into radia- 
tion the distribution of the particles in energy space will 
change. If monoenergetic electrons (8) are injected quasi- 
continuously and are confined in the radiating volume for 
a time long compared to the synchrotron loss timescale 
(~ 10 5 sec to allow cooling down to 7 ~ 1 with the pa- 
rameters considered here), the particle distribution will 
become 



N[p/) cx I7I 1 cx 7 



(19) 



This produces a radiation spectrum I v cx v 5 which is 
^ ' very flat. This is rather stable because the particles now 
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are isotropic in the comoving jet frame and do not have the 
large beaming factor as before but only the one from the 
bulk plasma motion, namely T. Together with the beam 
radiation one has a highly variable polarised v 03 compo- 
nent and a less variable or constant v~ - 5 power-law com- 
ponent with a different degree and angle of polarisation. 
This combination seems to fit to the observations quite 
well. A two-component model was already suggested by 
Qian ct al. (1991) to explain the drastic variations in the 
linear polarisation. 

Assuming that the beam is quasi-steady we also have 
to consider the inverse Compton scattering of the power- 
law component from the background jet by high energy 
beam particles. Since the secondary component radiates 
isotropically in the jet frame (also backwards), the beam 
electrons will have head-on collisions with those photons 
(radio to optical) . The photons scattered in the observer's 
direction will be boosted in frequency by a factor of order 
7o, thus producing a power-law radiation component from 
the optical, over the X-ray to the gamma-ray part of the 
spectrum. This emission is beamed towards an observer 
who also receives the radio to optical radiation, probably 
with a time lag. The timescale for the variability of this 
high frequency radiation is also determined by fluctua- 
tions in the beam structure and should therefore be given 
by eq. (14), i.e., it is of the same order of magnitude as 
for the variations in the radio-optical part. In fact this is 
seen in IDV-blazars (Wagner and Witzel 1995). Since the 
variability time-scale depends on the high Lorentz factor 
of the beam electrons and not on the bulk Lorentz factor 
of the background jet plasma, the constraints on the com- 
pactness of the emitting region are much less stringent. 

5. Summary and discussion 

We consider a model for BL Lac objects and quasars that 
show variabiliy on timescales less than a day. The basic in- 
gredient is a monoenergetic relativistic electron beam. We 
discuss the formation and stability of such an anisotropic 
distribution and find that it is stable as long as the mag- 
netic field along the jet is strong enough. The beam can 
propagate along the field lines without being scattered. 
When the phase velocity of the Alfven waves starts to de- 
part significantly from the velocity of light the beam will 
excite magneto-hydrodynamic waves. The particles will be 
scattered in pitch angle and start to emit synchrotron ra- 
diation at radio to optical frequencies. The electrons in the 
beam are then coupled to the background plasma of the 
jet. The spectrum of the monoenergetic particles which 
arc injected quasi-continuously evolves to a power-law. 
This gives rise to a second component in the radiation 
spectrum. Concerning the high brightness temperatures 
derived from variability arguments in IDV blazars, the 
upper limit on the compactness of the emitting region is 
relaxed by the Lorentz factor of the beam, which can be 
much larger than the bulk Lorentz factor of the jet. 



It seems even possible, when the Lorentz factor of the 
beam particles is of the order of 70 ~ 10 6 , that also highly 
variable emission in the TeV range observed in several 
AGN (e.g. Mkn 421; Gaidos et al. 1996) can be produced 
this way. There the observed timescales go down to 0.5 
hours. This is beyond the capabilities of models involving 
bulk relativistic beaming of intrinsically isotropic distri- 
butions of particles and radiation with particle accelera- 
tion at shocks. When the TeV emission comes from an 
anisotropic beam distribution the intrinsic length is only 
restricted to /' ^ 10 20 (At/l hr)(7 /10 6 ) cm. This reduces 
the opacity of 7 — 7 interactions below its critical value. 

For the objects considered here the Lorentz factor of 
the monoenergetic electron beam is in the range 10 2 — 10 4 . 
This is the basic input parameter. The model presented 
in this paper can account for several features observed in 
IDV blazars: 

1. ) the high variability brightness temperatures by rela- 
tivistic effects 

2. ) synchrotron radiation from the monoenergetic beams 
has an inverted spectral index as observed in the most 
variable sources 

3. ) the dominant timescale is quasi-stable, which can be 
understood in terms of the distance from the black hole 
where the standing shock developes 

4. ) together with the radiation from the reprocessed beam, 
a two-component model can be constructed, to explain 
drastic changes in the polarisation and its angle 

5. ) in such a combination I v oc azA 3 + bv~ o b one also 
expects the radio variations to be more pronounced at 
higher frequencies, which is seen in the modulation index 

6. ) since the synchrotron radiation of the beam reaches 
from the radio up to the optical, correlated variability 
(with the same timescale) is obviously explained 

7. ) a larger variability amplitude at optical frequencies 
compared to the radio is expected (and observed) 

8. ) correlations of the radio spectral-index with optical 
flux (flatter spectrum when optical flux rises) without any 
measurable time lag can be understood in terms of the 
varying v 1 / 3 emission, which gives the simultaneous fluc- 
tuations at radio and optical frequencies and also a spec- 
tral flattening (increase of spectral index) of the summed 
components 

9. ) because the degree of polarisation of synchrotron ra- 
diation from a monoenergetic distribution depends on the 
frequency (contrary to a power-law distribution) this also 
gives a key to the observed variations and the frequency 
dependence of the polarisation and its angle 

10. ) if the emission of the two components are polarised 
(nearly) orthogonal to each other, this goes along with 
anticorrelated variations of the total and the polarised flux 
(this is the case for BL Lac objects, where the magnetic 
field is perpendicular to the jet) 

11. ) if the polarisation of the two components is parallel, 
correlated variability of the total and polarised flux is ex- 
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pected (this is the case for quasars, where the magnetic 
field is oriented parallel to the jet) 

12. ) time lags between total and polarised emission may 
be attributed to the reprocession of beam particles into 
the power-law distribution 

13. ) variability at X-ray and gamma-ray wavelengths 
(with timescales comparable to the radio/optical) can be 
accounted for in head-on collisions of photons radiated 
backwards from the power-law component (comoving with 
the background jet) off beam electrons by inverse Comp- 
ton scattering. 

This model seems to be compatible with many prop- 
erties of IDV blazars. It is now necessary to work out the 
details of the acceleration and radiation processes and to 
compare them with specific objects. 
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